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The nonlinear I-V characteristics of mesoscopic samples contain parts which are linear in the 
magnetic field and quadratic in the electric field. These contributions to the current are entirely 
due to the electron-electron interaction and consequently they are proportional to the electron- 
electron interaction constant. We present detailed calculations of the magnitude of the effect as a 
function of the temperature, and the direction of the magnetic field. We show that in the case of a 
magnetic field oriented parallel to the sample, the effect exists entirely due to spin-orbit scattering. 
The temperature dependence of the magnitude of the effect has an oscillating character with a 
characteristic period on the order of the temperature itself. We also clarify in this article the nature 
of the electron-electron interaction constant which determines the magnitude of the effect. 

INTRODUCTION 

According to Onsager, the linear conductance G(H) of a conductor measured by the two-probe method must be 
an even function of the magnetic field H Q] : 

G(H) - G(-H) (1) 

This is a consequence of general principles: the time reversal symmetry and the positive sign of the entropy production. 
Therefore it holds in all nonmagnetic conductors. In a single particle approximation and at zero temperature the 
validity of Eq. ^ also can be verified using the Landauer formula for the conductance of a sample 

G= X^ |T4 ''' (H)|2 (2) 

and requirement of time reversal symmetry Ty(H) = T* (— H). Here Ty- is a scattering matrix between electronic 
channels labelled by the indices i and j (see for review of the subject Ref. Q). To verify the validity of Eq. 2] in 
mesoscopic samples one has to prove that ((G(H) — G(— H)) 2 ) = 0, which requires the cancellation of the odd-in-H 
parts of the Cooperon and Diffuson type of diagrams shown in Fig. lh,i,j. Here the brackets () denote averaging 
over random realizations of the impurity potential (we use standard diagram technique for averaging over random 
scattering potential |^). 

On the other hand there are no general principles preventing the existence of odd-in-H terms in the nonlinear I-V 
characteristics of conductors. In this article we study the quadratic-in-voltage V part of the I-V characteristics which 
can be represented as 

I {n l ) =V 2 [F (H) + F e (U)] (3) 

where F (H) and F e (H) are odd and even functions of H respectively. Since H is an axial vector and the current 
density, j, is a polar one, the function F (H) can be non-zero only in non-centro-symmetric media. It is important to 
■ study F (H) because of the fact that in the approximation of noninteracting electrons, F (H) = 0. It is particularly 
simple to verify this fact using the Landauer formula. Indeed, in the absence of the electron-electron or electron- 
phonon interactions, the total current through the sample will be the sum of contributions from different electron 
energies. Each of these contributions is an even function of H, and hence the total current will also be an even 
function of H. Thus the effect is entirely due to electron-electron or electron-phonon interaction. In contrast, the 
even-in- H function F e (H) can be described even in a single particle approximation (see, for example, Ref. where 
the calculations were done in the case of mesoscopic samples). 

In the case of pure bulk non-centro-symmetric crystals, and at high temperatures the effect described by F (H) 
has been investigated both theoretically and experimentally (see, for example, Ref. Q). In the case of chiral carbon 
nanotubes a classical theory of this effect was discussed in Ref. |8| • At high temperatures there are two contributions 
to the effect: 
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a) The first contribution is purely classical and it can be described in the framework of the Boltzmann kinetic 
equation: An electric field accelerates electrons creating a non-equilibrium distribution function. This non-equilibrium 
distribution function consists of two parts: an anisotropic- in-momentum part that is proportional to V and a quadratic- 
in-V contribution which is isotropic in the momentum. In isotropic media, the relaxation of this non-equilibrium 
distribution due to inelastic scattering processes yields no net current; however, in non-centrosymetric media and in 
the presence of the magnetic field, the inelastic relaxation rate has odd in electron momentum components, which 
give rise to the odd-in-H part of Eq. described by F G (H). 

b) The second contribution yj is due to a shift in the center of mass of a wave-packet during collisions. The 
description of these processes is beyond the classical Boltzmann kinetic equation. In non-equilibrium and non-centro- 
symmetric media and in the presence of the magnetic field these shifts take place in a particular direction determined 
by a crystal symmetry, and lead to an odd in H contribution in the net current through the sample. This contribution 
is similar to that discussed in the framework of the anomalous Hall effect 01- The above two contributions to F (H) 
are proportional to the inelastic electron relaxation rate l/r e , and consequently, they vanish at T = (there are of 
course contributions from the above effect to the I-V characteristics which at T = are proportional to higher powers 
of V). 

It has been pointed out in Refs. [ill fl2| that in mesoscopic metallic samples, where all spatial symmetries are 
broken, there is an odd-in-H contribution to Eq. [3] which survives in the limit T = and which therefore determines 
the magnitude of the effect at small T. This effect has been observed experimentally 0,0, ElEJEJ- As usual for 
mesoscopic effects, this contribution is due to random electronic interference. Therefore it exhibits random sample 
specific oscillations as a function of the external magnetic field, temperature and the electron chemical potential. The 
characteristic feature of the effect is that it is proportional to the amplitude of the electron-electron interaction rather 
than the scattering rate. The qualitative explanation of the effect is the following The linear in V mesoscopic 
fluctuations of the current density are due to random interference of electron waves travelling along different diffusive 
paths. Though the total current through the sample should be an even function of H, the local current densities 
contain a part which is odd-in-H. By the same token, there is a part of the electron density An(r, V, H) which is 
proportional to V and odd-in-H . In the presence of the spin-orbit scattering the applied voltage also induces local 
fluctuations of the spin density AS(r, V,H). These nonequilibrium densities create an additional random potential 
due to electron-electron interaction 

Au(r,V,H)=/3g ) / An(r,V,H) (4) 

and an additional exchange magnetic field 

h(r ! y,H)=/3g ) / AS (5) 

(1 2) 

were P^ff are interaction constants. We can then calculate a change of a linear conductance of a sample, AG = 
G({u(r) + Au(r), h(r)}) — G({u}), induced by a change of the scattering potential given by Eq. which gives us 

I (nl) =AG[V,n,{Au(r,V,H), h(r,V,H)}]V (6) 

Expanding Eq. EJwith respect to A«(r) , and h(r) we arrive at Eq. |j3 

In this paper, we present calculations of the magnitude of the effect as a function of the temperature and the 
magnetic field. We show that in the case of parallel magnetic field, the effect exists entirely due to spin-orbit 
scattering. The temperature dependence of the magnitude of the effect (both in parallel and in the perpendicular 
magnetic field) randomly oscillates with a characteristic period on the order of the temperature itself. We also take 
into account the exchange contribution to the effect and clarify the nature of the electron-electron interaction constant 
which determines the value of F (H.) in EqOH 



DIAGRAMMATIC CALCULATIONS OF THE MAGNITUDE OF THE EFFECT 

Let us consider a sample shown in the insert of Fig. 1 with the characteristic size L which is much larger than the 
electron elastic mean free path I. In this limit (I n i) = 0. Therefore we shall characterize the magnitude of the current 
by the variance (/„>)■ Before averaging over random potential configurations, to first order in the electron-electron 
interaction constant, the value of I n i is given by the diagram shown in Fig. la. After averaging, the quantity (Z^j) 
is given by diagrams shown in Fig. Id. In these diagrams, the solid lines which correspond to the electron Green 
functions carry frequencies on the order of the temperature T, so we can neglect frequency dependence of interaction 
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propagators. Thus one can introduce an additional scattering scalar potential given by Eqs. [5] substitute it into 
Eq. and arrive at Eq. |3] 

Generally 0^ff an< ^ @eff are phenomenological parameters represented by thick wiggly lines in Fig. la, d. At high 
electron density, electrons are weakly interacting. In this limit the diagram in Fig. la is reduced to those shown in 
Fig. lb, c. These diagrams correspond to the Hartree and Fock contributions respectively. This means that p e tf can 
be calculated to first order in perturbation theory with respect to the electron-electron interaction V(r) 

1- 



0$ f = V(0)-^V(k F (n -■*>)) 



0% = -V{k F {n-n')) (7) 

where V(q) is a Fourier transform of a screened Coulomb interaction. The bar denotes the average over the angle of 
the unit vectors n, and n'. To verify this fact one has to show that diagrams Fig.le,f,g contains only combinations 

(V(0) - iV(Mn-n'))) and (i^(M n-n'))) ■ 

The system of equations 1 11516171 is a generalization of that in Ref. where only Hartree term was taken into 
account. Eq. is typical for many effects in mesoscopic conductors with interacting electrons (see for example 
[l?l l2o| ). Usually, however, electron-electron interaction effects give small corrections to the conductance of good 
conductors with G 3> e 2 /h. In our case, the magnitude of the effect is proportional to Peff- 

To get Eq. |3|one has to expand the expression for the conductance in Eq. [SI with respect to Au(r) and h(r). To 
do so it is convenient to expand the potential [2l| 

oo 

Au(r) = ^Uimir) (8) 

i 

and the effective magnetic field 

oo 

h a (r) = ^2 h a ,iTii(r) (9) 

i 

in a complete set of orthogonal eigenstates rii(r) of the diffusion equation (J drnf(r) = 1): 

d 2 

D-^-m(r) = E ini (r) (10) 
err 

where Ei are the eigenvalues of Eq^J| and i labels the eigenstates. We assume boundary conditions, which correspond 
to zero current through a closed boundary, and n.j(r) = at the open boundary. Generally speaking the electron 
density Ara(r) and, consequently Au(r) contain all spatial harmonics, and the problem is similar to the sensitivity of 
the sample conductance to a change in the scattering potential Ait(r) considered in [22Tl23| . Thus we have 

In the absence of the spin-orbit scattering the second term in Eg II II is zero while dG/dui is an even function of 
H. Thus in this case the effect originates from the odd in H part of An(r, V, H). In the presence of the spin-orbit 
scattering the second term is nonzero and there is another contribution to the effect which comes from the odd-in-H 
part of h(r, V, H). Using Eq. II II and neglecting small correlations between AG, Au(r) and h(r) we get 



i fin An 

) ((fti, a (H) - ^, a (-H))(/i^(H) - hj,/9(-H))) (12) 



4V 2 ^ X du t duj 

dG dG 



— t 

AVI 



4^/2 dh la dhj 

i,j,a,0 ' J " 

Eq. H2l onlv contains correlation functions which can be estimated in a single particle approximation (in zero order in 
V(r)). This can be done in a standard way (see for example [Til H^. Eif ) by calculating diagrams shown in Fig.le,f,g. 
These diagrams contain ladder parts, shown in Fig.lij, which depend on the electron spin indices. After summation 
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over the spin indices, the diagrams shown in Fig.le,f,g only contain the blocks Pj(r, r 1; ei 2 , H,u) s ) (i — 1, 2, 3) described 
by the equations 

[-D (V + if A) 2 + 2 ^ + iw s si 5 nu; + r" 1 ) P a (r, r',w, H,w s ) + ^ (P x - P 2 ) = 5 (r - r') 

(-£> (V + if A) 2 + iw - iuJsSignuj + r" 1 ) P 2 (r, r',w, H,w s ) + ^ (P 2 - P x ) = (13) 
(-D (V + if A) 2 + iuj + J- + iu s signuj + t" 1 ) P 3 (r, r',w, H,cj s ) = 5 (r - r') 

Here t,;, is the spin-orbit mean free time. These equations account for the dependence of Pi on spin-orbit scattering 
and magnetic field. The magnetic field enters through the vector potential A and the Zeeman splitting lu s = g/isH, 
where /is is the Bohr magneton. The boundary condition for these diffusion poles at the insulating boundary is 
n (V + if A) Pi(r, r',ui, H,cj s ) = 0. In the case of ideal leads, when the electron diffusion coefficient in the leads 
is infinite, the boundary condition at the leads is Pi(r, r',oj, H,cj s ) = 0. The quantities Pj(r, r 1; ei2, H,w s ) enter 
expressions for /„; only in the combination 

J2 l=1 2 (P(r, r-^2,0, 0)Pi(ri, r,e 21 ,0, 0) - P(r, y iM2 , 2H,^)P(n, r,e 2 i, 2H,w fl )) + , . 

+P 3 (r, r^fi, w s )P 3 (r l5 r,e 21 ,0, u a ) - P 3 (r, r^eia, 3H,0)P 3 (ri, r,e 21 , 3H,0) 1 J 

It is convenient to choose the gauge nA = at insulating boundary. In this gauge, when calculating magnetic field 
dependence we can use standard perturbation theory with respect to the magnetic field. 

The sum over i in Eq. H2l converges quickly and, consequently, the main contribution comes from the zero-harmonic 
with i ~ (this fact is related to the long range character of the correlation function of the part of the electron 
densities which are proportional to V |l8| ). Consequently the approximation where only this zero harmonic of the 
potential 

tio(V,H) = f An(r,V,H)dr (15) 

V J v 

is taken into account gives a result valid in order of magnitude. In this formula, v is the volume of the sample, which 
reduces to A, the area of the sample, in the two-dimensional case. In the case when r s -C 1 ^> fif^ we also can 
neglect the exchange field h(r, V, H) . 

In the rest of the article we consider the case of a two-dimensional sample which has a thickness much smaller that 
its lateral dimension, L L z . In this case the results are different for cases of a perpendicular and a parallel magnetic 
field. 



Effect in a Perpendicular Magnetic Field 

In the case where the magnetic field is oriented perpendicular to the film, spin-orbit scattering can be neglected as 
long as t so ^> t. In this case the result depends on the relations between the sample size L, the magnetic length Lh , 
the dephasing length L v , and the thermal coherence length of normal metal Lt- It also depend on the nature of the 
leads to the sample. 

Let us start with the case of ideal leads when the electron diffusion coefficient in the leads is infinite, Dl — oo. 
Then at small temperatures Lt, L v >Lwe get 

tf m )-B*tfl,r^(*)\»)' (i6) 

where T = D/L 2 , D = vpl/2 is the diffusion coefficient inside a two-dimensional sample, B is a numerical factor of 
order one, and A is the area of the sample. 

The linear ii-dependence of /„; given by Eq. 1 161 holds at small magnetic fields when HL 2 <C $0j where $o = fie/c 
is the flux quanta. In the opposite limit HL 2 3> $o the function P D (H) exhibit random oscillations as a function 
of H with a characteristic period of order AH ~ <J>o /L 2 . To verify this one can show that the correlation function 
(F (H + AH)F (H)) decays at AH ~ $ /L 2 . Qualitatively the dependence P G (H) in the case Dl = oo is shown in 
Fig. 2a. 

The temperature dependence of P (H, T) is nontrivial because it exhibit random oscillations as a function of T 
which are shown qualitatively in Fig. 2b. To verify this one has to calculate two correlation functions 



(P o 2 (T))~(F o (0)F o (T)) 



(17) 
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and notice that they have the same temperature dependencies. It also follows from Eq. El that at T > D/L 2 : the 
period of oscillations is on the order of the temperature itself. Let us turn now to the case when the leads have the same 
diffusion coefficient as the sample, Dl = D. In this case H and T dependencies of F(H,T) are qualitatively shown 
in Figs. 2c, d. The main difference with the case Dr, = is that in this case F (T,H) exhibits random oscillations 
as functions of H and T even in the case Lh S> L and Lt 3> T. These oscillations are related to the existence of 
diffusive electron trajectories which leave the sample, travel in the leads and then come back. Their contribution to 
the total current is small, but their sensitivity to changes of H and T are so big that the derivatives 8F /dT — ► oo and 
8F /dH — + oo at zero temperature and magnetic field H — ► 0. Qualitatively the T and H dependencies of F a (T,H) 
are similar of those of the linear conductance 1221, 12J, |2a . 



Effect in a Parallel Magnetic Field 

When deriving Eq. EH we neglected the effect of Zeemann splitting of the electron spectrum because it gives a 
small contribution to the effect. However, when a thin enough sample is oriented parallel to a magnetic field H||, the 
orbital contribution of the magnetic filed is absent, and Zeemann splitting causes the dominant contribution to /„,;. 
However, in the absence of spin-orbit scattering, <$n(H||) , and consequently J n ; are even functions of H||. Thus in 
this case the effect is determined by the spin-orbit scattering rate 1/t so . In the presence of spin-orbit scattering the 
dashed lines in Fig.l correspond to the following expression 



1 



2'kvt 



Safins* 



1 



8itvt so 



(18) 



In the case of a weak magnetic field, /ib ■ H|| -C 1/t so ,T, and at T = 0, we get an expression for the odd-in-H|| 
part of the zero mode of the density fluctuation 



(i£(H,,)>=C - 



h 2 A 2 T 2 



0&) 



2 1 



1 



r(r + -3- 



i 



i 



2i^ + 2r — + 



/ fi B H 



V r 



(19) 



where C ~ 1. In the case of weak spin-orbit scattering T ^> r^ 1 we have (F 2 ) — T~ 7 . It is interesting that 
the amplitude of the effect is proportional to the amplitude of the spin-orbit scattering rather than it's rate. This 
expression holds as at HbHu < T. In the opposite limit we have (F 2 ) = T~ 5 which is independent of H||. In this 
region .F (H||) exhibit random sample specific oscillations as a function of Hu with a characteristic period of order 

T/fig. To get this result one has to calculate the correlation function ([F (H|i) — ^(Hii + AH||)] 2 ) and to see that 
it decays at fxHn > T. 

In the limit of a strong spin-orbit scattering T, /isHii <C 1/t so we have (F 2 ) ~ 1/t so , and the amplitude of the 
effect which decreases as r so increases. In the limit of strong magnetic field /isH|| > 1 we have (F 2 ) — t 2 /T 5 . In 
this case the period of the oscillations of F (H||) as a function of H|| is of order 1/t so . 



CONCLUSION 



We have shown that odd-in-H and quadratic-in-T^ part of the current is proportional to the electron-electron 
interaction constant. Thus detailed measurement of this effect should yield information about the strength of the 
electron interaction. On the other hand, in mesoscopic samples the amplitude of the current given by Eq. 3 is 
proportional to a random sample-specific sign. Thus it is unclear at the moment whether it is feasible to extract the 
sign of the electron interaction from measurements of the nonlinear current Eq. |3J 

At some level the effect considered above is similar to the effect of the interactional corrections to average con- 
ductivity of disordered metals These corrections originate from a correlation between random electron diffusive 
trajectories and the Friedel oscillations of the electron density in the presence of random potential. The difference is 
the following. In the case of Ref. electrons scatter on equilibrium Friedel oscillations of the electron density, which 
are even in H. At G ^ e 2 /h this leads to non-analytic but small and even in H corrections to the Drude conductance. 
The origin of our effect is the electron scattering on non-equilibrium fluctuations of the electron density. In this case 
the electron interaction determines the magnitude of the effect. The effect considered above is also different from 
classical effects the amplitude of which is proportional to the inelastic relaxation rate and which does not exhibit 
oscillations as a function of temperature, magnetic field and the chemical potential. 
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The magnitude of the effect discussed above decays as the temperature increases. If the crystalline structure of the 
material is non-centro-symmetric, at high enough temperatures, the T-dependence of F (H,T) is determined by the 
"classical" effects considered in Ref. U0@- However, if the structure of the pure crystal is centro-symmetric, then 
the only source of the effect at high temperatures is the non-centro-symmetric distribution of the scattering potential 
u(r). In this case the effect is of a classical nature. Namely, one should consider the classical motion of interac ting 
electrons in the presence of a frozen random potential, similar to what has been done for average quantities in |27|. 
This problem however, is beyond the scope of this article. 

Finally we would like to discuss a difference between our approach and the approach in Ref. 0]. At T = and 
in the absence of the electron interaction the Landauer formula in a combination with the random matrix theory has 
been a useful tool for describing the linear conductance of mesoscopic samples. In this approximation it can be derived 
from the Kubo formula [28j. However, even in this case the derivation can be carried out only in the case of ideal 
leads, Dl — oo, when incident and transmitted waves through the sample are well defined. At finite temperature, and 
in the absence of inelastic processes L <C equation still can be applied. However, to describe the temperature 
dependence of the conductance G(T) one has to know delicate properties of the energy dependencies of the matrix 
elements Ty(e). At even higher temperatures, L ^> y/Dr e , when inelastic scattering processes are significant, Eq. 
121 cannot be applied. One of the reasons for this is that the electron channels (and even their number) in Eq. [21 
are not well defined in this case. For example, the Landauer formula cannot reproduce the well-known Bloch T 5 
and T 2 temperature dependencies of resistivity of bulk metals associated with electron-electron and electron-phonon 
scattering [29j . It also can not reproduce electron interaction corrections to the conductance considered in jlflj . 

Sometimes the Landauer approach gives correct results for linear conductance of samples even in the case when the 
leads are not ideal and Dl < oo. Consider for example a constriction between two semi-infinite 3D metals and assume 
that the diffusion coefficient is independent of the coordinates. The Landauer formula still gives a correct result for 
the conductance of the sample in a single particle approximation because the conductance is determined by the part 
of the sample near the constriction. On the other hand, the magnetic field G(H + 8H) — G(H) and the temperature 
G(T + ST) — G(T) dependencies of the conductance in this case are determined by the interference of diffusive paths 
travelling on distances of order of the magnetic length Lsh and the coherence length Lst of the normal metal. At 
small AT and <5H these lengths are much bigger than the constriction size, and they diverge as 5H, ST — > 0. As a 
result, in the case of non-ideal leads, Dl < oo, and in the absence of inelastic phase breaking processes, the periods 
of oscillations of G(H,T) as a function of ST and SH decrease at small H and T, and the derivatives dG/dH and 
dG/dT diverge at T — ► and H — > 18]. To cut off these divergencies one has to take into account inelastic electron 
scattering processes. Such effects arc beyond the Landauer formula. 

The situation with non-linear parts of the I-V characteristics is more complicated. In the single particle approxima- 
tion the I-V characteristics can be expressed in terms of the e-dependence of the matrix elements Tij(e), where e is the 
electron energy. This procedure gives a result equivalent to that obtained in |5| using Keldysh diagram technique. In 
this approximation, however, F a {H) — 0. Generally speaking in the presence of the electron interaction the Landauer 
approach can not be justified even in the case of ideal leads, and even at T = 0. One of the reasons is that in the 
presence of electron-electron and electron-phonon inelastic processes allowed at V ^ 0, the electron channels are not 
well defined. In other words the voltage plays the role similar to the temperature, and at V ^ there exist effects 
which are similar to T-dependent interactional corrections to conductivity [19j. On the other hand, at small voltages 
these processes lead to the value of I n i proportional to a power of V greater than or equal to three. Then the question 
arises whether the Landauer scheme can be modified to describe the quadratic-in-V current of Eq. The authors 
of Ref. [TH ] introduced a concept of non-equilibrium capacitance which relates V and the total charge induced in the 
sample. On a phenomenological level this approach is similar to those presented in Ref. and in this article. We 
would however like to mention, the differences. The approach of Ref. [Ill ] corresponds to accounting for only the zero 
harmonics of Au(r) in Eq. Though this approximation is not exact, in diffusive samples and at T = it gives the 
correct order of magnitude of the effect. In ballistic quantum dots the mistake is much bigger. More importantly the 
approach of Refs. [ll| is restricted to the Hartree approximation and can not account for the exchange interaction. 
We would like to mention that before averaging over realizations of the scattering potential, the exchange terms in 
Eqs. 0] and [SI can change even the sign of the effect. 

In conclusion we would like to mention that there must also exist currents through the sample which are proportional 
to ff(VT) 2 and HVVT . 
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FIG. 1: Diagrams describing I„i, and (1^;). Thin thin solid lines correspond to the electron green functions, thick solid lines 
correspond to electron Green functions averaged over random realizations of the scattering potential, dashed lines correspond 
to the correlation function of the scattering potential given by Eq. 1181 . The fat wiggly line is the effective electron interaction 
interaction, thick wiggly lines correspond to the elecron-electron interaction V(r). Diagrams b,c corresponds to Hartree (b ) 
and Fock (c) terms of the perturbation theory with respect to V(r). 
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FIG. 2: Qualitative pictures of H- and T dependencies of Fo(H,T). Figs, a) and b) correspond to the case of ideal leads 
(Dl = oo), while Figs, c and d correspond to the case of diffusive leads Dl < oo . 



